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Abstract 

We give a gauge-invariant definition of the vortex surface in SU (N) Yang- 
Mills theory without using the gauge fixing procedure. In this construction, 
gauge-invariant magnetic monopoles with fractional magnetic charges emerge 
in the boundary of the non-oriented vortex surface such that the asymptotic 
string tension reproduces the correct iV-ality dependence. We show that gauge- 
invariant magnetic monopoles and vortices are simultaneously responsible for 
quark confinement in four dimensional spacetime based on the Wilson criterion. 
These results are extracted from a non-Abelian Stokes theorem derived in the 
previous paper. 
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1 Introduction 



In understanding the non-perturbative phenomena in the infrared sector of Yang- 
Mills theory [1] and QCD such as quark confinement [2], chiral symmetry breaking 
and U(1)a problem, some of topological configurations are believed to play the key role 
as the dominant dynamical degrees of freedom. Examples are magnetic monopoles, 
center vortices, instantons, merons, etc. Among them, chiral symmetry breaking and 
the Ua(X) problem can be explained e.g., by the Yang-Mills instantons 0,11], although 
magnetic monopoles are not excluded as their mechanisms, see e.g. [5]. 

Quark confinement is believed to be explained by the condensation of Abelian 
magnetic monopoles [6] and/or center vortices [HIE], since they realize the dual super- 
conductor picture of QCD vacuum as the most promising scenario for confinement. 
Magnetic monopoles are topological objects of codimension 3: points for D = 3 and 
closed loops for D = 4 in D-dimensional spacetime. On the other hand, center vor- 
tices are topological objects of codimension 2: closed loops for D = 3 and closed 
surfaces for D = 4 [HHHHIII]- For example, the vortex in D = 3 is in the first step 
identified with a closed thin tube of magnetic flux which can be thought of as the 
magnetic field generated by a toroidal solenoid in the limit of vanishing cross section, 
although the vortex must have a finite transverse extension to correctly reproduce the 
adjoint string tension and yield the finite action of vortices [12]. Recent lattice sim- 
ulations exhibit both Abelian magnetic monopole dominance p3| and center vortex 
dominance [9] for the string tension. Moreover, if either Abelian magnetic monopoles 
or center vortices in Yang-Mills theory are removed from the ensemble of configura- 
tions, confinement is found to be lost. Moreover, chiral symmetry breaking is also 
lost. It is known according to numerical simulations that magnetic monopoles and 
vortices are strongly correlated. Incidentally, merons [16J can be also a candidate for 
confiners jl] and have something to do with magnetic monopoles due to numerical 
simulations on a lattice [T7]. For recent reviews, see e.g., [H] for Abelian magnetic 
monopole and [15] for center vortices. 

Yang-Mills instantons and merons are respectively self-dual Euclidean and non- 
self-dual Euclidean/Minkowski solutions of the gauge covariant Yang-Mills field equa- 
tions derived from the Yang-Mills action. However, Abelian magnetic monopoles [18] 
or center vortices [7] in SU(N) Yang-Mills theory have been obtained as gauge fix- 
ing defects by a partial gauge fixing G = SU(N) — > H where the gauge degrees of 
freedom is used to transform the gauge field variable (link variable on a lattice) as 
close as possible to a subgroup H which is left unbroken: the maximal torus (Car- 
tan) subgroup H = [/(l)^" 1 for the maximal Abelian gauge [T8l[T9] or the center 
subgroup H = Z(N) for the maximal center gauge. Therefore, the current method 
of constructing Abelian magnetic monopoles and center vortices could not escape the 
charge of the gauge artifact. 

In a series of recent papers, we have succeeded to give a gauge-invariant description 
of the dual superconductivity in Yang- Mills theory in the continuum [2"3"}l2"4"j and on a 
lattice [231111123123122] by developing the approach founded in [2H2I1I22] • Especially, 
the Wilson loop operator is expressed exactly in terms of a gauge-invariant magnetic 
current [33] and the magnetic monopole can be defined in a gauge-invariant way 
according to a non- Abelian Stokes theorem for the Wilson loop operator [30jl31U32[l33] . 
These results enable us to clarify the role of magnetic monopole in confinement. 

In this paper, we give a gauge-invariant definition of a vortex which can play the 
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same role as the center vortex to sweep away distrust of gauge artifact. Indeed, this 
is achieved without relying on the (partial) gauge fixing such as the maximal Abelian 
gauge and the maximal center gauge. Consequently, the gauge- invariant magnetic 
monopoles emerge in the boundary of the closed vortex surface in D = 4. We show 
that both magnetic monopoles and vortices are necessary from a viewpoint of quark 
confinement based on the Wilson criterion such that the asymptotic string tension 
reproduces the correct iV-ality dependence. The result of this paper shows that the 
vortex and the magnetic monopole can be the alternative view of the one and the same 
fundamental dynamics in Yang-Mills theory defined in a gauge- invariant manner. 



2 Wilson loop operator and magnetic monopoles 

For the Yang-Mills connection (one-form) = £/^(x)dx^ = stf^ (x)T A dx M for a gauge 
group G, the Wilson loop operator Wc\s^\ along a closed loop C is defined by 



WcW] := tr 



exp \ ig 



/tr(l), (2.1) 



where 8? denotes the path-ordering prescription. The non- Abelian Stokes theorem 
(NAST) enables us to rewrite the Wilson loop operator into the surface integral 
form over the surface S bounding C (<9E = C). A version of the non- Abelian Stokes 
theorem without any path or surface ordering is known as the Diakonov-Petrov version 
[3"U] of a non- Abelian Stokes theorem. The Diakonov-Petrov version of NAST was 
originally derived in [30] for G = 577(2) case and later developed and extended to 
G = SU(N) case in [3U[32]. Moreover, it has been shown [33] that the Wilson 
loop operator is rewritten in terms of two gauge-invariant conserved currents, the 
"magnetic- monopole current" k and the "electric current" j, defined by applying the 
exterior derivative d, the coderivative (adjoint derivative) S and Hodge star operation 
* to /: 

k:=5*f = *df, j:=Sf, (2.2) 

where / is the gauge-invariant two-form defined from the SU(N) gauge connection 
<(x) = tf{x)T A by 

U„(x) =dn2ti(n(x)j*u(x)) - d„2tr(n{x)^{x)) 

+ 2tr( 2(iV iv ~ 1) ^- 1 n(x)[^n(x),^n(x)]), (2.3) 

using the su(N) Lie-algebra valued field n(x) = ua{x)Ta called the (normalized) 
color field defined by 



n(x) = n A (x)T A := J J^Lfe)^ ( x ), £(x) <= G/H, (2.4) 
with the stability group H specified later and H given by 

H:=AH = jZA 3 H v (2.5) 
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where Hj (j — 1, , ■ ■ ■ , r) are the generators from the Cartan subalgebra of su(N) 
(r = N — 1 is the rank of the gauge group G = SU(N)) and r-dimensional vector Aj 
(j = 1, , ■ • • , r) is the highest weight of the representation in which the Wilson loop is 
considered. Indeed, both currents are conserved in the sense that 5 k = and 5j = 0. 

For the Wilson loop operator in the fundamental representation of G = SU(N), 
H = U(N — 1) [33] and the Wilson loop operator is rewritten as 



WcW\ = /[tMOkexp 



<N- 1 
2N 



9 I f 

Jt, 



= /[rfMO] S exp|z^^^(A ; ,H E ) + z^^^( J ,iV s )|, (2.6) 

where is the product of the Haar measure on SU(N)/U(N — 1) over 

S: 

[dtit)h ■= II Mt(?))> (2-7) 

and S E is the (.D — 3)-form and iV E is the one-form defined by 

S s : = *de s A~ 1 = 5* OsA" 1 , iV E := <5e s A-\ (2.8) 

with the .D-dimensional Laplacian (or the d'Alembertian in the Minkowski spacetime) 
A := d5 + 5d and the two-form E called the vorticity tensor as an antisymmetric 
tensor of rank two: 

Q%{x) : = jf dS^{X{a))5 D {x - X(a)). (2.9) 

Note that k and H E are (D — 3)-forms, while j and iVs are one-forms for any D in D 
dimensional spacetime: 



(*,Se) := * / d D xA^ 1 -^-»(x)S^- to -»(x), := / <Px,f (x^x). 

(2.10) 



For G = SU(2), in particular, arbitrary representation is characterized by an 
integer or a half- integer J = |,l,|,2,|,---. The Wilson loop operator in the repre- 
sentation J of SU(2) obey the non-Abelian Stokes theorem: 

WcW\ = J [dfi(Oh exp {Ug(k, S E ) + ^(j, iV E )} . (2.11) 

This agrees with (12.61) for a fundamental representation J = ~ of SU(2). 
We focus on the magnetic contribution W™ defined by 



3 



For D = 3, it has been shown [33] that the magnetic charge q m defined by q m = J d 3 xk° 
for SU(2) obeys the quantization condition: 

q m = 4ag- 1 n, n G Z = {■■• ,— 2, — 1,0,+1, +2, ■ • ■ }. (2.13) 

This follows from the condition that the non-Abelian Stokes theorem should not 
depend on the surface E chosen for bounding the loop C, since the original Wilson 
loop is defined for the specified closed loop C. For an ensemble of point-like magnetic 
charges located at x = z a (a = 1, • ■ ■ ,n) 

n 

k°(x) = J2Ct {3) (x-Za), q a m = ^g~ 1 n a , n a EZ, (2.14) 

a=l 

we have a geometric representation: 

W™ = exp f; OM*a)j = exp |z| £ n a fi s (^ a )| , n a G Z, (2.15) 

where fis (x) is the the solid angle under which the surface E shows up to an observer 
at the point x. Therefore, a magnetic monopole with a unit magnetic charge q m = 
An g^ 1 in the neighborhood of the surface E gives a non-trivial factor exp[±z7r] = — 1 
to the Wilson loop operator W™, since ^ls{.z a ) = ±2n when z a is just below and 
above the surface E. This is a nice feature of magnetic monopole for explaining 
quark confinement based on the Wilson loop. 

For D = 4, however, we fall in a trouble about the magnetic monopole just defined. 
For D = 4, an ensemble of magnetic currents on closed loops C' a (a = 1, ■ ■ ■ , n): 

n - 

fc"(aO=£&/ dy^\x-y a ), q a m = ^g~ x n a , (2.16) 

a=l JC 'o. 

leads to [33] 

Wg = exp Uig q a m L(C' a , E)| = exp S^ni £ n a L(C' al E)| , n a G Z, (2.17) 

where L(C, E) is the linking number between the curve C and the surface E: 

L(C",E) = L(E,C") := / dy^r)E^y(r)), (2.18) 

where the curve C is identified with the trajectory of a magnetic monopole and the 
surface E with the world sheet of a hadron (meson) string representing a quark- 
antiquark pair. For D = 3 case, see Fig. HJ However, such magnetic loops carrying 
the magnetic charge obeying the quantization condition (12.131) do not give non-trivial 
contributions to the Wilson loop, since n a and L are integers. If the quantization 
condition (12.131) is true, the magnetic monopole can not be the topological defects 
responsible for quark confinement. In the following, we discuss how this dilemma is 
resolved. 
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I(C,S) Iff ,C) 

Figure 1: In D = 3 dimensional space-time, a non-trivial contribution to the Wilson 
loop operator comes from (Center panel) the linking L(C, C) between the closed loop 
C and a closed (magnetic current) loop C as a boundary of the vortex sheet S. (Left 
panel) intersection I(C = <9E, S) between the closed loop C and an open vortex sheet 
S bounding the closed loop C, or (Right panel) intersection J(E, C = dS) between 
the open surface £ bounding the closed loop C and a closed loop C . Three are 
equivalent descriptions: L(C, C) = I(C = dE, S) = /(£, C = OS). 



3 Magnetic monopole and vortex 

In the following, we extensively use the techniques developed by Engelhardt and 
Reinhardt [TT] in constructing a continuum analogue of the maximal center gauge 
and center projection, but from a different angle in this paper. 

We first consider the D = 4 case. Suppose that the magnetic current has the 
support on the closed loop C in D = 4 dimensions: 

kHx) = k"(x; C) := $ / dy^ix - y), (3.1) 

Jc 

where $ is a real number representing the magnetic flux carried by the magnetic 
charge q m to be discussed later in detail. The magnetic charge q m is defined by 

q m = J d 3 a^, (3.2) 

where x^ denotes a parameterization of the 3-dimensional volume V and d 3 a^ is the 
dual of the 3-dimensional volume element <fV 717273 : 

doc^^~ doc^^~ dx^^~ 

d 3 a^ := ^e At7l7273 rfV 717273 , d 3 cr 717273 := e Pl p 2 p 3 ^— ^ ^ da 1 der 2 d(T 3 . (3.3) 

First of all, we look for the field strength *f fiu (x; S) with the support S, a two- 
dimensional surface bounding the closed loop C, dS = C, so that (12. 2p holds: 

d v *r{x;S)=W{x;C'). (3.4) 

Such a solution is given by 

*f^( X ] 5) = $ J sds d 2 a^S 4 (x - x(a)). (3.5) 
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In fact, it satisfies the equation: 

d u *f^(x;S) =$ f d 2 a^d x j\x-x{a)) 
J s 

=$ f d 2 a^dl5\x-x{a)) 
•J s 

=<£>/ dx>"8 A {x-x{a)), (3.6) 

JdS=C 

where we have used the Stokes theorem in the last step. 

Next, we proceed to obtain the gauge potential b^(x; V) giving the field strength 
f^ u (x; S) just obtained, i.e., 

f^(x; S) =0M^ V) - d u b,{x- V), (3.7) 

such that b^x; V) has the support only on the open set V, the three-dimensional 
volume whose boundary is S: dV = S. Note that *f (J,v (x; S) is cast into 

*f^(x; S) =$ / d*a^ K dl5\x - x{a)) 

JV:dV=S 

=<$>e^±e Pm l v . dv= /° lll2l3d ^ x ~ 

=e^d* a [<£>±6p m [ _ d 3 a'^6\x-x(a))} 

= / d?~a?5\x - x{a))\ - (a <- 0)}, (3.8) 

2 Jv.dV=S 

where we have used the Gauss (Stokes) theorem in the first equality. Therefore, the 
gauge potential b^x; V) is determined up to a gauge transformation: 



b^x; V) = $ / d 3 a^5 4 (x - x(a)). (3.9) 



V:dV=S 



This corresponds to an explicit (singular) gauge field representation of an ideal vortex 
configuration in 4 space-time dimensions [H] . 

In D dimensions, the magnetic current k is (D — 3)-form with the support on the 
(D — 3) -dimensional subspace C' D _ 3 in D-dimensional spacetime: 



k{x) = k(x; C' D ^) := $ 6 d D ~ 3 x(a)5 D (x - x(a)). (3.10) 

Then it is not difficult to show that an ideal vortex configuration should read 

b^x; V) =$ f dP^a^ 5 D (x - x(a)). (3.11) 

JVD-i:dVD-l=SD-2 

The gauge field b^x; V) for the vortex is not unique as mentioned above. Actually, 
the ideal vortex field b^x; V) can be gauge transformed to a thin vortex field a^(x; S) 
which has the support only on the boundary Sd-2 = 9Vn-i of Vd-i- 

b^x- V) = a^x; S) + iU(x; V)d^(x; V), (3.12) 



1 The precise position of the open set V is irrelevant for the value of the Wilson loop as shown 
below. In fact, the open set V can be deformed arbitrarily by singular gauge transformations in such 
a way that its boundary dV representing the position of the magnetic flux of the vortex is fixed. 
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so that a^(x] S) carries the same flux located on S as that carried by b^(x; V). We 
adopt the gauge transformation: 

U(x; V) = exp [i$£l v (x)) , (3.13) 

where Qy(x) is the solid angle taken up by the volume V viewed from x: 

w^L^m^- A Sr (3 - 14) 

with Ad_\ being the area of unit sphere S D ~ l in D dimensions. The solid angle 
Qv{ x ) defined in this way is normalized to unity for a point x inside the volume V. 
Note that the solid angle is defined with a sign depending on the orientation of V as 
rays emanating from x pierce V. A deformation of V keeping its boundary S fixed 
leaves the solid angle invariant unless x crosses V. When x crosses V, the solid angle 
changes by an integer. H 

Now we show that the thin vortex field is obtained in the form: 

a^x- S) = $t d D - 2 ~a^ d x D(x - x{a)), (3.15) 

J Sd-2— dVjj-i 

where D(x — x(a)) is the Green function of the .D-dimensional Laplacian defined by 

-d^D{x - x{a)) = 5 D (x - x(a)), (3.16) 

with the explicit form: 

D( X -y) = - D ,T /2 -j m . (3.17) 

V yJ 4% D / 2 {(X - y)2p-2)/2 V ! 

The ideal vortex field b^x; V) has support only on V and hence it vanishes outside 
the vortex V, i.e., b^(x; V) = for x V. Therefore, outside the vortex V, i.e., 
x G M D — V, the thin vortex field a M (x; S) is the pure gauge due to (13. 12ft : 

a^x- S) = -iU(x; V)d fJ U\x ] V), x £ V, (3.18) 

which implies the vanishing field strength f^ u (x) = outside the vortex V. This 
is reasonable, because the magnetic flux is contained in the vortex sheet S. The 
magnetic field computed from the curl of a^(x; S) can be localized on C = dS. 

The derivation of (I3.12p and (13. 15j) is as follows. By using the fact that the solid 
angle is rewritten as 

n v (x) := f d^a^Dix-xia)), (3.19) 

JVn-i 

2 Whether the flux of a vortex is electric or magnetic depends on the position of the D — 2 
dimensional vortex surface S in Z?-dimensional spacetime. For example, in D = 4 the vortex defined 
by the boundary of a purely spatial 3-dimensional volume V carries only electric flux, which is 
directed normal to the vortex surface S = dV. On the other hand, a vortex S = dV defined by a 
volume V evolving in time represents at a fixed time a closed loop and carries the magnetic flux, 
which is tangential to the vortex loop [TT] . 
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we obtain 

iU(x; V)d p U\x; V) = $d p n v (x) = $ f dP'^d^Dix - x{a)). (3.20) 

JVd-i 

We have the decomposition: 

dpdv = 5 pu d 2 - (5^d 2 - d^dy) = 5 pu d 2 - ^e ppa pe uaa pd p d a , (3.21) 

which is used to rewrite the pure gauge form into 
zU(x;V)d fl U^(x;V) 

=$ f d D ~ 1 a p d 2 D(x - x(a)) - $ f d D ~ 1 aJ-e ppa , 3 e uaa pd p d (T D(x - x(a)) 

JVo-i JVd-i a 

= -$/ d D - l an5 D (x-x(a)) / d D ~ 2 d pX d x D{x - x(a)), (3.22) 

JVd-i JdV D -\ 

where we have used the definition of the Green function and the Stokes theorem in 
the last equality. In other words, the thin vortex field a M (x; S) is the transverse part 
of the ideal vortex b p (x; V). 

Finally, we find that the surface integral of f(x; S) over E bounded by the Wilson 
loop C is equivalent to the line integral of a(x; S) along the closed loop C: 



f(x;S) = idb(x;V) 

b(x; V) 



<9£=C* 

S)+iU(x; V)dU\x]V)} 

c 

a(x;S), (3.23) 

ic 

since the contribution from the last term iU(x; V)dW(x; V) = Qd p Qy(x) vanishes 
for any closed loop C. Then the line integral is cast into 

j> dx^a^x; S) = J d 2 a ufl (x)d v a tl (x; S) 

=$ / d 2 a vp {x) f d D ~ 2 &^d u d x D(x - x(a)) 

=<&—/" d D ~ 2 a pa {x)e pvpa I d 2 (j a p(x)6 Pj \ a pd v dxD(x - x) 

4 JT JdVo-i 

=$zl I d D ~ 2 d a p(x) [ d 2 a a/3 (x)d 2 D(x-x) 

2 JT JdV D -! 

+ $ / d D ~ 2 apa(x) I d 2 a aP (x)d a d a D(x - x) 

J S J &Vd — i 

=$~ / d D - 2 b a(i {x) f d 2 a Q( s(x)6 D (x-x) J (3.24) 

2 JT JdVn-i 

where we have used the fact that the second term vanishes due to the Stokes theorem 
ddV = for obtaining the last result. Therefore, the line integral is rewritten as 

I dx»a p (x; S) = $/(£, S D - 2 = dV D - X ), (3.25) 
J c 
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I(C,V) KE,S) 



Figure 2: In D = 4 dimensional case, (Center panel) Linking L(C, S) between the 
closed loop C and a closed 2-dim. vortex surface S bounding a 3-dim. vortex V. 
(Left panel) Intersection I(C = <9£, V) between the closed loop C and the vortex V. 
(Right panel) Intersection I(E, S = dV) between a minimal surface £ bounded by 
the closed loop C and a closed vortex surface S in 4-dimensional space-time. They 
are equivalent: L(C, S) = I(C = <9E, V) = /(£, S = dV). 



in terms of the intersection number /(£, S) between the world sheet £ of the hadron 
string and the vortex sheet S. It is known that the intersection number I(E, S") is 
equal to the linking number L(C = <9£, S = dV) between the Wilson loop C and the 
vortex sheet S, see Fig. |2j 

/(£, 5 D _ 2 = 9Vb_0 =^(C = <9£, 5 D _ 2 = 9V D _!) 

=i / rf D " 2 a a/3 (x) / dV /3 (x)5 D (x-x). (3.26) 

2 JS 

Thus the thin vortex contributes to the SU (2) Wilson loop in the fundamental 
representation: 



cxp 



{ 2 9 L f( x ' s ) 



exp 



i^g*L{C,S = dV) 



z L ( c,s=dV)^ z , = jh*. (3.27) 



In general, z is a complex number of modulus one, i.e., an element of U(l). If the 
magnetic charge obeys the Dirac quantization condition: 

$ = 2>Kg- 1 n (n G Z), (3.28) 

then z reduces to the center element Z 2 of SU(2): 

z = e inn = ±1, zl E Z 2 . (3.29) 

The quantization condition (I3.28P will be called the fractional quantization condition 
which happens to agree with the Dirac one for SU(2), which is realized as a special 
case of the general quantization condition for SU(N) discussed later. If the magnetic 
charge obeyed the quantization condition 

$ = Aitg-hi (n 6 Z), (3.30) 

then z would be trivial, i.e., z = 1. Such a vortex can not give a non-trivial contribu- 
tion to the Wilson loop. Therefore, the thin vortex carrying the fractional magnetic 
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Figure 3: The schematic view of three-dimensional slice for D = 4 spacetime dimen- 
sions. A vortex surface S intersects another surface S' which includes £ bounding 
C. The vortex surface S is not globally oriented, but consists of patches of different 
orientation. Boundaries of such patches are tagged by magnetic monopole (indicated 
as a thick line around the equator of the surface S). 

charge yields a center element under the fractional magnetic charge quantization 
(I3.28p . For D = 4, V is a three-dimensional volume and its boundary S is a closed 
two-dimensional surface S = dV. If S is an oriented closed surface, then its boundary 
dS = C is empty, and hence the magnetic current k^{x\C) does not exist in this 
case, since its support C = dS = ddV vanishes. Therefore, for the non-vanishing 
magnetic current k^(x; C) to exist in the boundary C = dS of S, the vortex surface 
S = dV must be non-oriented. Is there any relationship between the non-orientedness 
of the vortex surface and the fractional magnetic charge (13.281) ? We will try to answer 
this question in the next section. 

4 Non-orientedness of the vortex surface and frac- 
tional magnetic charge quantization 

The above considerations suggest that in the continuum theory a smooth vortex 
surface S consists of surface patches S n of different orientations S = U n S n , that is to 
say, the vortex surface S is not globally oriented and that the magnetic monopole loops 
fc™ emerge at the boundaries dS n of vortex surface patches S n where the magnetic flux 
direction on dS n is defined by the orientation of the patch S n (and vice versa). Then 
we can understand only a fraction of the elementary magnetic charge dictated by the 
Dirac quantization condition is carried by the magnetic flux $ on an isolated patch 
S n of the vortex sheet S. Hence the magnetic current of proper magnetic charges 
satisfying the quantization condition (I3.30P is reproduced only when the open vortex 
surface patches S n are glued together to obtain a non-oriented closed vortex surface S 
by combining their magnetic loops k™ together. See Fig. [3] In fact, if two patches are 
glued together such that the surface orientation does not change across the boundary 
(the surface is globally oriented), the magnetic current at the boundary precisely 
cancel. For D=4, thus, the thin vortex defined on a non-oriented closed surface S if 
linked to the Wilson loop C = <9E, gives a non-trivial contribution to the Wilson loop 
average. See Fig. [5] 

Engelhardt and Reinhardt [TT] have shown that the continuum center vortex con- 
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figurations generate the Pontryagin index ±1/2 as self-intersections of the vortex net- 
work. It is well known in topology that the self-intersection number of closed, globally 
oriented two-dimensional surface in IR 4 vanishes. This implies that the Pontryagin 
index vanishes for globally oriented vortex surfaces. Conversely, non-orientedness 
of the surfaces is crucial for generating a non-vanishing topological winding number. 
Therefore, the global non-orientedness of the vortex surfaces is necessary to generate a 
non-vanishing Pontryagin index originating from the vortex configuration. Moreover, 
it is well known in topology that the number of intersection points of two closed two- 
dimensional surfaces in M 4 is even. This implies that the number of self-intersection 
points of a closed surface is even, because the self- intersect ion number is defined by 
simply intersecting the surface with another surface infinitesimally displaced from 
it, i.e. a framing of the surface [TT]. Thus, even if each self-intersection point of 
an center vortex surface configuration gives a contribution ±1/2 to the Pontryagin 
index, the number of such contributions is even and the Pontryagin index is integer- 
valued. Moreover, they have pointed out that a non-zero Pontryagin index requires 
the existence of magnetic current. In the lattice study [TO] , in fact, the orientability 
of the vortex surfaces was studied, with result that these surfaces are non-orientable 
and have non-trivial genus in the confinement phase. These considerations connect 
the gauge-invariant magnetic monopole and the meron by way of a vortex and they 
become topological objects simultaneously responsible for confinement. 

5 iV-ality dependence of asymptotic string tension 

Now we show that the fractional charge quantization is crucial to understand the N- 
ality dependence of asymptotic string tension. The asymptotic string tension depends 
only on the iV-ality of the quark charge (or the group representation of the Wilson 
loop). This is because particles of zero iV-ality can never bind to a particle of 
non-zero iV-ality to form a color singlet and hence they can never break the string 
connecting two non-zero iV-ality sources. 

Gluons belong to the adjoint representation of the gauge group and they have zero 
iV-ality, so that all center elements are mapped to the identity. Therefore, gluons can 
not break the string formed between a pair of quark and antiquark in the N and 
N* representations of SU(N). Only other quarks, or other particles in the non-zero 
iV-ality representations, can do that. The Wilson criterion for quark confinement 
should be understood to imply the linear potential rising indefinitely in the limit in 
which the masses of any matter particles of non-zero iV-ality are taken to infinity. 

For example, in SU (2) gauge theory, the center group is Z,2 and the representations 
can be divided into two classes: J = |, |, • ■ ■ (half-integer) with iV-ality one, and 
J = 1, 2, ■ • ■ (integer) with iV-ality zero. The asymptotic string tension of iV-ality=0 
( J = 1, 2, • ■ • ) must be zero, while the non-vanishing asymptotic string tension of all 
iV-ality=l ( J = 2)1)"'") m ust be the same: 




3 The iV-ality k of a given representation is obtained by the number of boxes in the corresponding 
Young diagram, mod N. If M(g) is the matrix representation of a group element g in a representation 
of iV-ality k, and z € is an element in the center, then M(zg) = exp(2irink/N)M(g). 
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This expected result is easily derived from the considerations given in this paper: The 
vortex contribution to the Wilson loop operator in the representation J reads 



exp 



exp[iJg<5>L(C,S = dV)} 



z 



L(C,S=dV) 



(5.2) 



where 



z :- 



JJg<Z> 




(5.3) 



provided that the magnetic flux or the magnetic charge obeys the same fractional 
quantization condition (13.281) . For 7V-ality=0, z is trivial and the string tension is zero. 
For A r -ality=l, zeZ 2 has the non-trivial value —1 when the vortex sheet S intersects 
the world sheet £ and the string tension is non-zero. Thus the above argument 
reproduces correctly the iV-ality dependence of the asymptotic string tension. The 
SU (3) case will be discussed later. 



6 Consistency with lattice magnetic charge 

We consider how the above argument is consistent with the results of numerical 
simulations on a lattice, e.g., the magnetic monopole dominance in the string tension 
|26j . For this purpose, it is instructive to recall the definition of the magnetic charge 
and the magnetic current on a lattice for U(l) gauge theory [HI]. The basic idea is: if 
a magnetic monopole is located inside an elementary spatial cube on the lattice, then 
the enclosed magnetic charge can be determined by measuring the total magnetic flux 
through the surface (a set of plaquettes) of this cube. For instance, the magnetic flux 
$jfc(x) through the surface P x ,jk of a plaquette with an area e 2 lying in the j/c-plane is 
related to the phase of the plaquette variable U x jk = exp[ige 2 Fjk(x)] = exp[ig&jk(x)] 
with Fjk being the component of the magnetic field in the direction perpendicular to 
the jfc-plane. 

The plaquette variable Up — exp[ig$p] is given by the product of the four oriented 
link variables Uf = exp[i9e] G (7(1) around the boundary of a plaquette P: gQp = 
Yji&op^e where —n < 9g < tt. Hence, the phase (angle) associated with a given 
plaquette variable Up — exp[ig& P ] must satisfy 

- 47r < g§p < Air. (6.1) 

However, it should be remarked that the plaquette variable Up = exp[ig<&p] G U(l) 
is a periodic function of g$p with a period 2ir. Therefore, the physical flux should 
be determined from gQp, modulo 27rn, since the plaquette value remains unchanged 
by shifting g$p by a multiple of 2ir. According to DeGrand and Toussaint [35], we 
decompose this angle into two parts: 

g&P = 2nn P + g$ P , n P = 0, ±1, ±2, —n < g® P < n, (6.2) 

which covers the whole range of $p (16. ip . If we add up the plaquette angles $p of the 
six plaquettes bounding an elementary cube c, we will obtain a vanishing result, since 
each link is common to two plaquettes which give rise to the sum of two phases of the 
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common link variable, equal in magnitude but of opposite sign, i.e., J2p<=dc^P = 0- 
Therefore, the magnetic charge Q m as the magnetic flux through the closed surface 
S = dc bounding the elementary cube c is given by 

2n 

Qm= E $P = - E — n P- (6-3) 

P£S=dc Pes=dc 9 

Thus the magnetic charge is a multiple of If a magnetic monopole is located in an 
elementary cube, then at least one of the plaquette angle must be larger in magnitude 
than it, so that there is a Dirac string (line) crossing the corresponding surface. The 
Dirac string can be moved around by making a large gauge transformation such that 
a particular link variable is mapped out of the principal value [— 7r, 7r]. But the net 
number of such strings leaving the elementary volume will not be affected by the 
gauge transformation. Moreover, the number of magnetic monopoles contained in a 
volume V is given by the sum of the magnetic monopole numbers of the elementary 
cubes making up the volume V. 

Next, we consider a lattice expression for the components of the magnetic current 
k p defined by k^ = \e pvpa d u F pa . For example, the first component of the current 
k := [k\,k 2 , k 3 ) is written in the form: 

h = d 2 F M + d 3 F 42 + «9 4 F 23 = (d 2 , da, <9 4 ) • (F 34 , F A2 , F 23 ). (6.4) 

Therefore, integrating k\ over the volume of an elementary cube with edges along the 
2,3 and 4 directions is equivalent to computing the flux of F := (F 34 , F 42 , F 23 ) through 
the surface of this cube due to the Gauss theorem: \ c dx 2 dx 3 dx^k\ = J s=9c dS ■ F = 
J2p£S=dc&P — —J2pes=dc^f n P- m order to compute the three components of k, we 
need to calculate the flux through the plaquettes P of three cubes having one link 
directed along the 4-axis. This computation is carried out in a completely analogous 
way as described in the above for the magnetic flux. Hence, by construction, each 
component of the magnetic current k^(x) on the lattice will be multiples of ^ when 
measured in lattice units: 

2-7T 1 2tx 
k i( x ) = — ^ W^^vO^) = — d v *n^ v (x), n pa (x) = 0, ±1, ±2, (6.5) 
9 * 9 

where d v is the (forward) lattice derivative d p f(x) := [f(x + eji) — f(x)]/e. The mag- 
netic current satisfies the topological conservation law d'^k^lx) = and constitutes 
the closed loop in the dual lattice where d' p is the backward lattice derivative defined 
by ff,J(x) := [f(x) - f(x - e£))]/e. 

Thus, it happens that the definition of the magnetic monopole due to DeGrand 
and Toussaint for the U(l) gauge theory is consistent with the observation made for 
the vortex in the continuum formulation for SU(2). The gauge-invariant magnetic 
charge defined on a lattice reduces to the above one due DeGrand and Toussaint 
by taking a special gauge. Hence the result of [26] is consistent with the observation 
made in this paper. 
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7 577(3) case 



In the case of SU(N) gauge group, it has been shown [33] that the magnetic charge 
q m measured by the Wilson loop is subject to the quantization condition: 



2tt 2N 



n, n G Z, 



(7-1] 



which is analogous to the Dirac type, but different from it. Suppose that the magnetic 
flux $ obeys the fractional quantization condition: 



, 2tt 2N 

Then its contribution to the Wilson loop is written as 



cxp 



where 



IN -I 



ig\ 



2N Jz 



f(x; S) 



cxp 



z 



L(c,s=av) 

N > 



(7.2) 



(7.3) 



zn := exp < i\ 



<N- 1 
2N 



-g$ 



D 2mf N n 



(7.4) 



The center of SU(N) is Z N . For z N l to belong to the center of SU(3), i.e., Z 3 3 
{1, e 27ri//3 l, e 47r4//3 l = e~ 2m / 3 l}, fa must take the fractional number except for a trivial 
case /jv = 0: 

/s=i, | (75) 
just as the aforementioned SU(2) case with the center subgroup Z 2 9 {1, e m l = —1}, 



/a 



(7.6) 



This may be related to the fact that each self-intersection point of a vortex surface 
contributes ±k/N (k = 1, ■ ■ ■ , N — 1) to the Pontryagin index 

Finally, we consider the relationship between center vortices and our gauge-invariant 
vortices for SU (N) case. We introduce TC given by [33] 



H := A H = ^AjHj, 



(7.7) 



where Hj (j 



r) are the generators from the Cartan subalgebra (r = N — 1 is 



the rank of the gauge group G = SU (N)) and r-dimensional vector Aj (j = 1, , • • • , r) 
is the highest weight of the representation in which the Wilson loop is considered. 
For SU(2), every representation is specified by a half integer J: 



Ai = J — — ,1,— ,2, ■•• , iJi = — -, 
2 2 2 



(7 
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and the fundamental representation J = | of ££7(2) leads to 

For SU(3), the highest weight of the representation characterized by the Dynkin 
indices [m, n] is given by 

, . / m m + 2n\ 

A = (A " A2) = U-^rj' (710) 

The explicit form of H, for the fundamental representations [m, n] reads using the 
diagonal set of the Gell-Mann matrices A3 and A§: 

1., „ . 1 /2m + n —m + n —m — 2n\ 
H = -(AxA 3 + A 2 A 8 ) = -diag (—3— , — 3 , 3 ) , (7.11) 

and we enumerate all fundamental representations 3: 

[m,n] = [l,0]:A=(i^=J :=!/!, H = idiag (| ^) , (7.12a) 

Kr>] = [-1,1]:A=(^,^):=, 2 , « = ^ (t> I t) ' (7 ' 12b) 
lm,n] = [0, -1] : A = (o, ^) := , 3 , H = idiag (^, =i, |) = (7.12c) 

and their conjugates 3*. Therefore, the center vortex is obtained by replacing the 
magnetic flux $ with the diagonal matrix 4-kH,: 



N - 1 

^^T^^^n, (7.i3) 

which indeed leads to the non-trivial elements of the center group Z^: 



IN -1 
z N =e W {i x l—— g ^ 

-fexp(47riW) = e 27TifN l G Z N , f N = k/N (k = 1, ■ ■ ■ , N - 1), (7.14) 

where /jv = corresponds to the trivial element 1. Thus the center vortices are 
replaced by our gauge-invariant vortices carrying the fractional magnetic flux /at$. 

If this observation is correct, the magnetic part of the SU(N) Wilson loop with 
an additional fractional factor f^: 



W™ = ex P {igf N ^^^(k,Zx)}, f N = k/N (k = l,--- ,N-1) (7.15) 

will reproduce the string tension for any k (k = 1, • • • , N — 1), just as confirmed 
for SU(2) gauge group in [26J. This is a conjecture derived in this paper. This is 
consistent with the center vortex mechanism for quark confinement. However, it will 
be rather difficult to identify the vortex structure in £77(3) case in numerical simula- 
tions. For the vortex surfaces for SU(N), N > 3, may branch and the superimposed 
magnetic fluxes in general also modify the type of vortex flux, i.e., its direction in 
color space. 
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8 Conclusion and discussion 



In this paper, we have given a gauge-invariant definition of the vortex surface S 
such that the gauge-invariant magnetic monopole emerge in the boundary of the 
surface S. We do not use any gauge fixing for obtaining the relevant vortex and 
magnetic monopole, such as maximal Abelian gauge and maximal center gauge which 
are extensively used in the conventional approaches. These results are obtained from 
the non- Abelian Stokes theorem for the Wilson loop operator [33] • 

For the topological defects such as magnetic monopoles or vortices to give a non- 
trivial contribution to the Wilson loop average, thereby yielding the area law, we have 
shown that the magnetic current must carry the fractional magnetic flux (charge) 
compared with that dictated by the ordinary quantization condition. Under this 
identification, our gauge-invariant vortices can play the same role as the center vor- 
tices in giving a non-trivial contribution to the Wilson loop average. Hence, the 
gauge-invariant vortices are expected to reproduce a number of nice results obtained 
so far by center vortices [To] . 

For center vortices, it was pointed out [UJ that the global non-orientedness of 
the center vortex surface is necessary to generate a non-vanishing Pontryagin index 
from the vortex configuration in D = 4 dimensional spacetime. The existence of the 
fractional magnetic flux is consistent with the non-orientedness of the closed vortex 
surface in D = 4 where the smooth closed vortex surface must consist of the surface 
patches of different orientations so that the magnetic flux in the boundary of an 
isolated patch inevitably becomes fractional. 

In this paper, however, we have considered the infinitely thin magnetic flux tube. 
Therefore, the asymptotic string tension is correctly understood to reproduce the 
correct iV-ality dependence of the string tension. In order to see the behavior of the 
string tension in an intermediate range of distance, e.g., Casimir scaling or Sine-law 
scaling of the fc-string tension, we need to consider the magnetic flux tube with a 
finite thickness at around 1 fermi and the associated thick vortex. This issue will be 
discussed in a subsequent paper. 

Finally, we consider what happens to vortices as the number of colors N is in- 
creased. It is known [36] that instanton effects (of order exp(— 1/g 2 )) vanish exponen- 
tially exp(— N) for large N, since g 2 is of order 1/N, or A := g 2 N is to be fixed in the 
large N expansion. Therefore, instanton gas disappears in the large N limit and pre- 
dictions that depend on thinking about instantons and an instanton gas (even a dense 
gas) will not be correct. The argument for quantization of the topological charge Qp 
of instanton in QCD starts with the boundary condition: the gauge field approaches a 
pure gauge at infinity, ^ — > Gd^G" 1 as \x\ — > oo. Thus, in a quark confining theory 
this boundary condition is not reasonable and the conclusion based on it are likely to 
be wrong, provided that QCD remains a confining theory as N — > oo, as argued by 
't Hooft [37J. Do vortices get suppressed like instantons? If so, how does one explain 
confinement using vortices for large N. We can argue that vortices are not suppressed 
even in the large N. According to [35] . quark confinement tells us that the vacuum 
cannot be regarded as being mostly pure gauge. In fact, the expectation value of the 
Wilson loop in any state that is mostly pure gauge will (just as in perturbation the- 
ory) not show a confining potential. For any approximation that yields confinement 
will have to include fields that fluctuate at infinity, corresponding to a vacuum that 
is not mainly pure gauge. As argued in the text, our vortex picture is accompanied 
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by merons with half-integral topological number Qp and the meron field is written as 
srfp, = which is not the pure gauge. Therefore, the existence of vortices are 

not in conflict with the large N expansion. This suggests that vortex gas will not be 
suppressed and survive in the large TV limit to be responsible for confinement. More 
details will be given elsewhere. 
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